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■^ We consider the type I see-saw model with two right-handed neutrinos and a nor- 

mal neutrino mass hierarchy and impose a zero coupling between the right-handed 
neutrino mainly responsible for the atmospheric neutrino mass and the electron 
neutrino. We derive a master formula which relates see-saw input parameters in a 

_G one to one correspondence with physical neutrino observables. Using the master 

formula we search for simple ratios of couplings consistent with current data on 

Qh neutrino mass and lepton mixing. We discover a minimal predictive example in 

which the right-handed neutrino mainly responsible for the atmospheric neutrino 

mass has couplings to (u e , u^, v T ) proportional to (0, 1, 1) and the right-handed 

*vj neutrino mainly responsible for the solar neutrino mass has couplings to (V e , v^, v T ) 

$> proportional to (1, 1, 3) or (1, 3, 1), with a relative phase rj = =Fvr/3, providing the 

link between leptogenesis and CP violation in neutrino oscillation experiments. 

We show how these patterns of couplings could arise from an A^ family symmetry 

^O model of leptons which predicts all the PMNS parameters in terms of the neutrino 

■^j- mass ratio 1112/1713, corresponding to Tri-bimaximal-Cabibbo mixing, accurate to 

one degree, with the prediction 5 ~ ±7r/2. 
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1 Introduction 

The type I see-saw mechanism [TJ provides a beautiful understanding of the smallness 
of neutrino masses as being due to the heavy right-handed Majorana neutrino masses. 
Furthermore it strongly suggests that light neutrinos are also Majorana. However, de- 
spite its attractive features, the see-saw mechanism provides no understanding of the 
observed bi-large lepton mixing. It also provides no insight into either the ordering (i.e. 
normal or inverted) or the mass scale of the neutrinos (i.e. the mass of the lightest neu- 
trino). Moreover the see-saw mechanism is difficult to test experimentally, unless the 
right-handed neutrino masses are at the TeV scale [2], and typically contains more pa- 
rameters than physical observables [3] . The larger number of see-saw parameters means 
that the neutrino Yukawa couplings cannot ever be determined from physical neutrino 
observables. Finally, apparently "unnatural" cancellations can take place when con- 
structing the effective neutrino mass matrix from the see-saw parameters. 

One attractive idea which avoids the last problem of "unnatural" cancellations in 
the see-saw mechanism is that the right-handed neutrinos contribute sequentially to the 
light effective neutrino mass matrix with hierarchical strength, leading to a normal mass 
hierarchy of physical neutrino masses m 3 3> m 2 ^> mi. This is consistent with recent 
Planck results which provide no evidence of quasi-degeneracy and lead to a bound jl] , 

J2 m i< 0.23eV (95%; Planck + WP + highL + BAO). (1) 

i 

The idea of such a sequential dominance (SD) [5J of right-handed neutrinos is that one 
right-handed neutrino of mass M atm dominantly contributes to the see-saw mechanism 
and is mainly responsible for the atmospheric neutrino mass ms, while a second sub- 
dominant right-handed neutrino of mass M so \ is mainly responsible for the solar neutrino 
mass 777.2, with a third almost decoupled right-handed neutrino of mass M dec being mainly 
responsible for the lightest neutrino mass vn\. Moreover, in the limit that the lightest 
neutrino mass mi — > 0, the almost decoupled right-handed neutrino becomes irrelevant 
and the see-saw mechanism reduces to the two right-handed neutrino case with a normal 
neutrino mass hierarchy, which reduces the number of see-saw parameters as we now 
discuss. 

The type I see-saw mechanism with effectively two right-handed neutrinos [B] has 
fewer parameters and predicts one massless neutrino, leading to simplified formulas for 
neutrinoless double beta decay and the sum of neutrino masses relevant for cosmology. 
Texture zeros can further reduce the number of see-saw parameters [HI El EH EH] • With 
one texture zero [7J [TU] , the 3x2 neutrino Yukawa matrix contains 7 real parameters, 
after charged lepton phase rotations, which is exactly the same number of parameters 
as non-zero physical neutrino observables, comprising 2 real positive masses m%, m$, the 
3 mixing angles #23, 613, 612, the oscillation phase S and the single Majorana phase j3. 



In the case of the two right-handed neutrino model with a normal hierarchy, if the 
"dominant" right-handed neutrino has zero coupling to the lepton doublet containing 
the electron, this was shown over a decade ago to bound the reactor angle 6*13 < m^lrnj, 
[10]. We shall refer to this zero coupling as a "dominant texture zero". For example, 
such a "dominant texture zero" could appear in the (1,1) entry of the neutrino Yukawa 
matrix, similar to the quark and charged lepton Yukawa matrices, where a (1,1) texture 
zero it is a key ingredient of the Gatto-Sartori-Tonin (GST) [12] and Georgi-Jarlskog 
(GJ) relations [T3]. Recently, Daya Bay [U] and RENO |l5j have measured #13 ~ 0.15 
near the upper limit of the bound 6*13 < m 2 /m 3 [3 [TD]. The fact that this bound is 
saturated could be indicative of a particular underlying pattern of Yukawa couplings, 
and this provides part of the motivation for revisiting the minimal see-saw model. In 
particular we shall derive a new "master formula" which relates physical observables to 
Yukawa couplings, providing the necessary tools for such questions to be addressed. 

The measurement of a non-zero reactor angle excludes tri-bimaximal (TB) mixing 
[IB"] . It is nevertheless convenient to express the solar, atmospheric and reactor angles 
in terms of (small) deviation parameters (s, a and r) from TB mixing [T7] (for a related 
parametrisation see [H 



1 1 r 

sin0i 2 = -Wl + s), sin0 2 3 = -Wl + a), sin6» 13 = — . (2) 

For example, following the Daya Bay and RENO results, it was pointed out that the 
lepton mixing angles are consistent with the so called Tri-bimaximal-Cabibbo (TBC) 
ansatz [19] where 



s = 0, a = 0, r = A, (3) 

corresponding to 

sin6>i2 = -/=, sm6 23 = —=, sin #13 = — =, (4) 

with A = 0.225 being the Wolfenstein parameter, yielding the TBC angles, 

012 = 35.26°, 023 = 45°, 013 = 9.15°. (5) 

The TBC lepton mixing matrix is given to second order by [19], 

Vtbc^ -^(i + Ae*) ^(l-|Ae i5 ) ^(1 - ^A 2 ) 

TBC mixing thus assumes maximal atmospheric mixing and trimaximal solar mixing. 
For the normal neutrino mass hierarchy, some of the global fits, reviewed in [20] , prefer 
an alternative ansatz which we refer to as TBC2, corresponding to, 

s = -A 2 , a = -A/2, r = A. (7) 



+ (9(A 3 ). (6) 



This TBC2 ansatz has the feature that the atmospheric mixing angle is in the first octant 
and the solar mixing angle is somewhat less than its tri-maximal value as preferred by 
some of the latest global fits at the one sigma level, while TBC mixing generally remains 
valid at the three sigma level or better for all the global fits [20]. We emphasise that 
the issue of whether the atmospheric angle is maximal or not remains unresolved, with 
SuperKamiokande having a mild preference for a non-maximal atmospheric angle (21] 
while T2K prefers a maximal atmospheric angle, although both preferences are less 
than one sigma and hence not statistically significant. This situation is expected to be 
clarified in the near future. Returning to our analysis, the ansatze we have used will 
serve as an approximate guide in searching for particular patterns of Yukawa couplings. 
Once a particular pattern is identified, it then takes a life of its own, independently of 
the TBC or TBC2 ansatz, and leads to predictions for lepton mixing independently of 
either of these two ansatze which were used to discover the pattern in the first place. 

The present paper divides into two parts. In the first part we shall perform a 
model independent analysis of the two right-handed neutrino model with a normal 
neutrino mass hierarchy and a single "dominant texture zero" in the electron component 
of the Yukawa matrix [7J [10]. This model provides an attractive and rather minimal 
extension of the Standard Model consistent with all current neutrino data, with a one 
to one correspondence between the Yukawa couplings and the physical observables. 
We shall derive a master formula comprising exact analytic results which express the 
Yukawa couplings in terms of general values of the physical masses and mixings. These 
results in principle do not rely on the SD assumption and are valid even for large 
values of reactor angle. However we show that, in practice, the observed reactor angle 
$13 ~ m 2 /m 3 implies and is implied by SD, where the "texture zero" in the electron 
component becomes the "dominant texture zero" . We also confirm that a second texture 
zero is not possible, at least without invoking charged lepton corrections which we do 
not consider in this paper. In fact the two right-handed neutrino model with two 
texture zeros [8j [9j [11] was recently shown to be not viable for the normal hierarchy case 
123. 24 1 . which provides a further reason why we mainly focus on the one texture zero 
case HI [ffil HE] here. 

The second part of the paper is concerned with applying the results obtained in the 
first part to model building. An important step in this direction is to use the master 
formula to derive exact analytic results for the ratios of Yukawa couplings in terms 
of physical neutrino observables, and to use these results, together with the current 
best fit values of neutrino masses and mixing angles, to search for simple patterns of 
Yukawa couplings. We find that the atmospheric neutrino mass can have couplings to 
(V e , z/^, v T ) proportional to, for example, (0, 1, 1) or (0, 1, 2), which are strongly dependent 
on the value of the CP violating oscillation phase 5. The right-handed neutrino mainly 
responsible for the solar neutrino mass can have couplings to (u e , u^, u T ) proportional to 
(1,4,2) or (1,1,3) or (1,3,1), which are strongly dependent on the value of both the 
oscillation phase and the Majorana phase and differ from the traditional constrained 



sequential dominance (CSD) values [25]. In particular we focus on CSD3 defined to be 
based on the atmospheric couplings (0, 1, 1) and the solar couplings being either (1, 1, 3) 
or (1,3,1), with a relative phase difference r\ = =F7r/3. We show how these patterns 
of couplings could arise from an A^ family symmetry model of leptons which predicts 
all the PMNS parameters in terms of the neutrino mass ratio 7712/1713, corresponding 
to approximate Tri-bimaximal-Cabibbo mixing to an accuracy of one degree with 5 ~ 
±7r/2. The model provides a link between leptogenesis and CP violation in neutrino 
oscillation experiments. 

We remark that, although the two right-handed neutrino model with one texture 
zero case has been discussed before [7J EH E], such studies were performed about a 
decade ago when the mixing angles were not known so precisely, so it is timely to revisit 
them, in the light of the first measurement of the reactor angle by Daya Bay and RENO 
which approximately saturates the "dominant texture zero" bound #13 < m 2 /7n 3 [71 [TU]. 
Moreover, as indicated above, we present new and exact analytic results which directly 
relate Yukawa couplings to physical neutrino observables. Previously in [7J \10\ only 
approximate relations were given, while in [TTJ exact relations between Yukawa couplings 
and physical neutrino observables were only given parametrically via the complex angle 
z arising from an orthogonal parametrisation. Here we give a simple derivation of exact 
analytic relations between Yukawa couplings and physical neutrino observables via a 
master formula which does not rely on any parameterisation. We also emphasise a 
further important difference between the present paper and the earlier ones is that we 
know now that the bound #13 < m^JTrto, is approximately saturated, and in this paper we 
seek a simple answer to this question in terms of ratios of Yukawa couplings with simple 
phase relations. The fact that we find a simple answer, namely CSD3, is very gratifying, 
since it opens up new possibilities for model building as well as linking leptogenesis to 
the PMNS phases. We go on to construct an explicit A4 model of leptons for one of the 
examples that completely determines the PMNS mixing matrix including all the phases. 

The layout of the remainder of this paper is as follows. In the first part of the paper 
we perform a model independent analysis. In section[2]we briefly show how the two right- 
handed neutrino model emerges as the limiting case of a three right-handed neutrino 
model with sequential dominance. In section [3] we discuss the two right-handed neutrino 
with a normal neutrino mass hierarchy and a texture zero in the electron component 
of either one of the two right-handed neutrinos, and derive a master formula relating 
see-saw parameters to physical parameters for this case. Using the master formula we 
show that SD emerges as a consequence of bi-large mixing and the observed reactor 
angle, with the texture zero being associated with the dominant right-handed neutrino. 
In section [4], using the master formula, we calculate ratios of Yukawa couplings and 
see-saw phases, as functions of the PMNS phases, for particular mixing patterns such 
as TBC and TBC2 which are consistent with present data. 

In the second part of the paper we discuss the implications of the results for model 
building. After a brief review in section [5] of the status of indirect models, in section [6] 



we scrutinise the results obtained in the previous section, looking for simple possibilities 
which could form the basis of new indirect models. We discuss three such examples in 
some detail, checking their viability using exact numerical results for PMNS parameters 
and leading order approximate analytic results in section [7j We also discuss the link 
between leptogenesis and PMNS phases in section[8J In section[9]we then focus on CSD3 
based on the alignment (0, 1, 1) for the atmospheric column and (1, 1, 3) or (1, 3, 1) for 
the solar column with a relative phase difference rj = =F7r/3 and show how these follow 
from Ai family symmetry model of leptons based on F-term vacuum alignment and 
spontaneous CP violation. The model predicts all the PMNS mixing parameters in 
terms of the neutrino mass ratio, corresponding to approximate Tri-bimaximal-Cabibbo 



mixing to an accuracy of one degree with 5 ~ ±7r/2. Section 10 concludes the paper. 



2 Two Right- Handed Neutrino Model as a Limiting 
Case of Sequential Dominance 

The starting point for our analysis is the see-saw mechanism in the flavour basis where 
the charged lepton mass matrix Me is diagonal with real positive eigenvalues m e , m M , m T 
and the three right-handed neutrino Majorana mass matrix M R is also diagonal, with 
real positive eigenvalues, M atm , M so \, Md ec , 



M E = m M , M R = M sol . (8) 





We write the neutrino Dirac mass matrix as 

m e,atm m e ,sol l m e,decJ \ 
l» D = I m £atm ^sol [^decl = ( ™ftm ™s^l [™dcc] ) > ( 9 ) 

m r,atm m -r,sol [ m -r,dccJ / 

in the convention where the effective Lagrangian after electroweak symmetry breaking, 
with the Higgs vacuum expectation value (vev) inserted, is given by 

C = -E~ L M E E R - VZm D N R - l -W R M R N R + H.c. , (10) 

where vl = (v e , v^, u T ) are the three left-handed neutrino fields which appear together 
with E L = {e L ,fi L , t l ) in the lepton doublets L = (L e , L M , L T ) and N R = {N atm , N soh N dec ) 
are the three right-handed neutrinos and we have defined the three Dirac column vectors 

aS 7n atmi m sol? m dec- 

The term for the light neutrino masses in the effective Lagrangian (after electroweak 
symmetry breaking), resulting from integrating out the massive right handed neutrinos 



i.e. the see-saw mechanism with the light effective neutrino Majorana mass matrix 



in' = m D M R ~w 
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Sequential dominance (SD) then corresponds to the third term being negligible, the 
second term subdominant and the first term dominant: 



l atm 7 



l D T 



M t 



> 



m£,m D 



T 



sol "'sol 



atm 



M 



sol 



> 



m£„7?.> 



D i 
dec "'doc 



M, 



dec 



which immediately predicts a normal neutrino mass hierarchy, 
m 3 > m 2 [ > mi ] , 



(12) 



(13) 



which is the main prediction of SD. 

We have labelled the dominant right-handed neutrino and Yukawa couplings mainly 
responsible for the atmospheric neutrino mass m^ as "atm", the subdominant ones 
mainly responsible for the solar neutrino mass m-2 as "sol", and the almost decoupled 
(sub-sub-dominant) ones mainly responsible for m\ as "dec". Note that the mass or- 
dering of right-handed neutrinos is not yet specified. We shall order the right-handed 
neutrino masses as Mi < M 2 < M 3 , and subsequently identify M atm , M so \, Md ec with 
Mi, M2, M3 in all possible ways. 

It is clear that in the limit that mi — > then the sub-sub-dominant right-handed 
neutrino and its associated couplings labelled by "dec" decouple completely and the 
above model reduces to a two right-handed neutrino model. In that limit we simply 
drop the third terms [in square brackets] in EqsJ8 13 in anticipation of this. 



3 The Two Right- Handed Neutrino Model with Nor- 
mal Hierarchy and Dominant Texture Zero 

3.1 Derivation of the Master Formula 

Without assuming SD, we write the see-saw matrices in a simple notation as, 



m 



D 




M, 



R 



Y 
X 



(14) 



where we have written the complex Dirac masses as a, b, c, d, e, / with d = and the 
real positive right-handed neutrino masses as Y, X. We are in a basis where the charged 



6 



lepton mass matrix M E is diagonal with real positive eigenvalues ordered as m e , m^, m T . 
Later on we shall show that the observed reactor angle 6^3 < rri2/m 3 implies and is im- 
plied by SD, where the "texture zero" in the electron component becomes the "dominant 
texture zero", i.e. we will later identify d = rn^ Sitm = and Y = M atm in the notation 
of the previous section. However, to begin with, we do not presume this identification. 
We also allow the real positive right-handed neutrino masses Y and X to have either 
mass ordering, i.e. namely Y = Mi or Y = M 2 where Mi < M 2 . Hence the texture zero 
in EqjM] covers both cases with a texture zero in the electron component, 





m v = x x , xx. (15) 



We now derive a master formula in which the see-saw parameters in Eq|l4] may be 
related to the low energy physical observables such as neutrino masses, mixing angles 
and physical phases. Integrating out the right-handed neutrinos the see-saw formula 
gives the approximate light effective Majorana neutrino mass matrix, 

m v = ~m D M^ 1 m DT . (16) 

in the convention where the effective Lagrangian is given by q 

£ = -E~ L M E E R - l -vZm v v c L + H.c. . (17) 

Performing the transformation from the flavour basis to the real positive mass basis by, 
V El M e V Er = diag(m e , m M , m T ), V„ L m u Vj L = diag(mi, m 2 , m 3 ), (18) 

the PMNS matrix is given by 

U = V El V} l . (19) 

Since, as before, we are in the basis where the charged lepton mass matrix M E is already 
diagonal, then in general V El can only be a diagonal matrix, 



(20) 



consisting of arbitrary phases, where an identical phase rotation on the right-handed 
charged leptons V Er = P E leaves the diagonal charged lepton masses in M E unchanged. 







/ g*</>e 








Ve l 


= Pe= 1 





g*^M 









V 





g#r 



1 Note that this convention for the light effective Majorana neutrino mass matrix m v differs by an 
overall complex conjugation compared to some other conventions in the literature. 



In practice the phases in Pg are chosen to absorb three phases from the unitary matrix 
VJ and to put U in a standard convention, 



U = VP 
where, analogous to the CKM matrix, 



C12C13 



•S12C13 



S13C 



V 



-S12C23 - C 12 S 13 S 2 3e tS C12C23 - Si 2 Si 3 S23e l<5 C13S23 

S12S23 - C 12 S 13 C23e l5 -C12S23 - Si2Si 3 C 2 3e 1 ' 5 C13C23 



(21) 



(22) 



and the Majorana phase matrix factor is, 



/ 



P 



■Pi 
e 1 ^ 

■f>2 

e l 2 



\ 



(23) 



V I) 

From Eqs p|T9pO| we find, 

WPeTT^PeU* = diag(mi, m 2 ,m 3 ). 
Then using Eqj2l] we find 

V^ PETn v ' PeV* = Pdiag(mi,m2,m 3 )P = diag(0, e i/3 m2,m 3 ), 



(24) 



(25) 



where we have dropped m\ = (since we are considering the two right-handed neutrino 
model with a normal mass hierarchy) and defined (3 = fa- Using the see-saw formula, 
Eqfj~6| this may be written, 



-P E m u M^m u P E = 1/diag(0,e^m 2 ,m 3 )^ 



or, 



m D 'M- R x m D ' A 



where, 



,-D' 



m = IPeVH 



D 



Vdiag(0,e^m 2 ,m 3 )V T , 



J4e 

e^" 

e^ 



a 
e b 
f c 



(26) 



(27) 



(28) 



Finally, expanding Eqj27] using Eqs. 14 and 28, we arrive at our "master formula" 

a 2 ab ac 

ab e 2 + b 2 ef + bc ) = e^m 2 V a2 V 02 + m 3 V a3 Vp3, (29) 

a/3 
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ac ef + bc f 2 + c 2 



where we have defined 

a = — ==-, b= — =^, c= — ^, e = — =-, f= — ==-. (30) 

Vx y/x Vx Vy Vy v; 



The "master formula" in Eq.29 is a very useful equation, since it relates the re- 
phased Dirac mass matrix parameters a, b, c, e, f, scaled by the positive square roots of 
the right-handed neutrino masses X, Y, which are the inputs of the see-saw mechanism 



see Eq. 14), to physical low energy observables. The combinations of see-saw input 



parameters appear on the LHS of Eq 29 while the physical neutrino masses, mixing 



angles and phases appear on the RHS of Eq. 29 



In terms of parameter counting, we see that the RHS of the "master formula" Eq.29 



with V in Eq.22, involves 7 non-zero physical parameters, comprising 2 real positive 
masses m 2 ,m 3 , the 3 mixing angles #23, #13, #12, the oscillation phase 5 and the single 
Majorana phase (3, which determine exactly the 5 complex parameters a, b, c, e, / on 
the LHS which therefore involve only two independent physical phases. We shall refer 
to a, b, c, e, / as the "physical see-saw parameter combinations" . Note that the 5 input 
complex Dirac masses a, b, c, e, / also only depend on two physical phase differences 
which are left invariant under charged lepton re-phasing as discussed further later. 

3.2 Sequential Dominance and the Dominant Texture Zero 

So far we have considered the two right-handed neutrino model with a single texture zero 
in the electron component, assuming a normal neutrino mass hierarchy, and have derived 



a "master formula" in Eq.29, relating see-saw parameters to physical observables. As 
a first application of this formula, we show that SD is a consequence of the assumed 
texture zero d = which turns out to be the "dominant texture zero" . 
Writing out the "master formula" EqJ29J explicitly for all (a, (3) gives, 

a 2 = e ifi m 2 V^ + m 3 V^, (31) 

~ab = e i/s m 2 V 12 V 22 + m 3 V 13 V 23 , (32) 

ad = e l/5 m 2 V 12 V 32 + m 3 V 13 V 33 , (33) 

ef + bc = e ip m 2 V 22 V 32 + m 3 V 23 V 33) (34) 

e 2 + b 2 = e il3 m 2 V 2 \ + m 3 V 2 2 3 , (35) 

P + c 2 = e^m 2 V 3 2 2 + m 3 V 3 2 3 . (36) 



We can solve for the individual parameters in terms of the above combinations, 
(abf 



a 2 : 
(ad) 2 



(37) 



c 2 = ^-, (31 

~e 2 = (e 2 + b 2 )-b 2 , (39) 

P = (f + £ 2 )-5 2 . (40) 

Using Eqs.31 36 and TBC mixing in Eqj6j we estimate, to leading order in A, 



a 2 ~ e*^, (41) 

~ab ~ e^ + e-^, (42) 

~ac ~ -e^ + e"^, (43) 

e~ 2 + & 2 ~ ^, (44) 

f + c 2 ~ f . (45) 



Since the mass ratio may be approximated as m2/m^ m 6\z ~ \/y/2, from Eqs.41 45 
conclude that, 



we 



e + /r» {a + b + cY (46) 



since the LHS is of order WI3, while the RHS is of order mi. Using Eq.30, we see that 



Eq.46 just corresponds to the SD condition in Eq|l2"1 where we identify the see-saw 



matrices in Eqjl4l with the notation in Section |2j 

n \ f m . n \ 

(47) 




Y \ ( M atm 



X ) \ M sol 

from which we see that the texture zero in the electron component is in fact the "dom- 
inant texture zero", since we identify mf atm = and Y = M atm . From our present 
perspective, where the reactor angle is measured, we see that SD and the "dominant 
texture zero" implies and is implied by the data on the mixing angles which is approx- 
imated here by TBC mixing. However, historically, a decade before the measurement 
of #13, SD was simply postulated, along with the "dominant texture zero". This was 
then shown to lead to the bound on the reactor angle #13 < m 2 /m 3 using approximate 
analytic formulas for the mixing angles and phases expressed in terms of input see-saw 
parameters pO, \10\ . 
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3.3 Excluding a second texture zero 



In this subsection we use the "master formula" in Eq, 29 to search for simple patterns of 
neutrino Dirac masses, or equivalently neutrino Yukawa couplings, which could motivate 
some future flavour model, for example based on a discrete family symmetry. To address 
such questions, it is relevant to take the ratios of the see-saw parameters in the same 
Dirac column, so that the right-handed neutrino masses cancel. From Eqs 31 36 and 

, V23V12 



Eqj30]we find, 
e 

b 

z 2 = - 
a 



z 3 



c 
a 



__ V13V22 

V33V12 — V13V32 
e l/3 m 2 V 12 V 22 + m 3 V 13 V 23 

e il3 m 2 V 2 2 + m 3 V 2 3 
e l/3 m 2 V 12 V 32 + m 3 Vi3V33 



(48) 
(49) 

e^m 2 V 2 2 + m 3 V^ (50) 

Note that Z\ is independent of neutrino masses (and phases). It is also worth emphasising 
that, using Eq.30 the moduli of these ratios are equal to the moduli of the original Dirac 
mass matrix elements, 



\ z i\ 



|e| 1 1 |fe| . . \c\ 

i/r ' w ' i°r 

For the case of TBC2 mixing in EqjTj we find, 
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Xe 
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V6 2 

~(1 - -A + -Xe* 5 + -A 2 - -AV 5 ) 
y/3 y 2 2 2 4 ' 

^(1 + -A--A 2 ). 
y/2 l 2 2 ; 



(51) 

(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 
(60) 
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Figure 1: The left and right panels show |z 2 | and |z 3 | (respectively) as functions of the physical phases 
8 and (3 evaluated using Eqs |49|50| and using the TBC2 ansatz in Eq|4j working to second order in A. 
From these plots we see that a second texture zero corresponding to cither |xa| = or \z^\ = is not 
possible. Similar results are obtained for TBC mixing. 



We shall approximate the mass ratio as: 
m 2 



= e 



m 3 



7i + 3 A 



0.176. 



(61) 



According to the global fits, which are compared in [20J, the mass ratio lies in the one 
sigma range m 2 /m 3 = 0.17 — 0.18. 

Using these approximations we perform an expansion of the ratios including terms 
of order A 2 , which we use for our numerical results. However, since the full expansion is 
not very instructive, we only write out the terms to order A below: 



e 



±[l-A(l + 2e- Jd )] + (9(A 



Z2 



a 



l + ^e 
\/2 



" (/3+<5) + ^X(B -C) + C(A 2 ) 



c 
a 



-1 + -=e 
V2 



-i(fi+s) + ^ X {B + C) + 0{\ 2 ) 



where we have introduced the complex notation, 



B 

C 



2e - i (/ 3 + 5 ) _ g e -i(2/3+3<5) 



3x/2e-^+ 25 ) + 



w 



-i(p+6) 



(62) 

(63) 
(64) 

(65) 
(66) 
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We can use these results to investigate the viability of a second texture zero. Ap- 
parently a second texture zero with \a\ = is inconsistent with Eqsj63p4} Also either 



e| = or \f\ = is inconsistent with Eq,62. A second texture zero with either |6| = 



or \c\ = would correspond to either \z%\ = or \z 3 \ = 0. This looks unlikely to be 
possible from the leading order results, but requires numerical investigation to be sure. 
We find that, for both TBC mixing and the more accurate TBC2 parameterisation given 
by Eqs J52|J60| to second order in A, either \z 2 \ = or \z 3 \ = are not possible. This is 
illustrated for the second case of TBC2 mixing in Figjl} In summary, we conclude that 
all the following two texture zero cases are not viable, 



m D 








(67) 



4 Searching for simple ratios of Yukawa couplings 
and phases 

We shall write the neutrino Dirac mass matrix, or equivalently the neutrino Yukawa 
matrix, in some arbitrary phase basis (but in the diagonal right-handed neutrino and 
charged lepton mass basis) as, 

(68) 

where the phase (fi a is not physical and can always be removed by charged lepton re- 
phasing, while the two physical phases which are left invariant under charged lepton 






ae* 


|e|e^ e 


\b\e i<t>b 


|/| e % 


| c |g#c 



re-phasing are given, using Eq. 30 , by 



r] 2 = 4>b - (fie = arg b - arg e, r] 3 = (fi c - (fi f = arg c - arg /. (69) 



Note that the columns of Eq.68 may be interchanged since the right-handed neutrino 
masses are arbitrary, but the texture zero is always in the electron component of the 
"dominant" column as shown previously. In order to determine the two physical phases 
r)2 and rh we first need to determine the second order variables on the left hand side of 
Eqs 31 36 From these we can then determine all the quadratic variables, a 2 ,b 2 ,c 2 ,e 2 ,f 2 . 



For example b 2 = (ab) 2 /a 2 . There is a sign ambiguity in taking the square root of these 
quadratic variables in order to find a,b,c,e,f which we need to discuss. We first note 
that the see-saw mechanism is invariant under a simultaneous change of sign of a, b, c, 
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or (independently) e, /, which corresponds to an independent undetermined sign in the 
each of the two columns of the Dirac mass matrix. Without loss of generality we shall 
use this sign freedom in the Dirac columns to fix b and e to be, 



arg 6 



-argfr 2 , 



arge 



arge 2 , 



(70) 



so that these phases lie in the restricted range — 7r/2 to tt/2, rather than the full range 
— 7r to 7T over which argfr 2 and arge 2 vary. Then the phase of c is determined using 



formulae which include the information about relative phases of b and c in Eqs,31 



36 



1 _ 2 ( ab 

arg c = - arg c — arg - 

2 \ ac 



(71) 



(where the relative phase of b and c is meaningful since these variables are defined such 
that the charged lepton phases are fixed). Having determined the phases of e, b and 
c, we can then determine the phase of / using a formula which takes into account the 
relative phase of / in the term involving ef + be in Eqs,31 36 which, after squaring and 
rearranging, yields the formula, 



arg f = — arg e — arg b — arg c + arg - (ef + bc) — e f —be 



'.2~2 



(72) 



Having thus correctly determined the phases of b,c,e,f, then 772, 773 may be determined 
using Eqj69] together with Eqs |70|7T|72 



Since we have excluded the possibility of a second texture zero in the previous sub- 
section, here we now consider other types of simple Yukawa coupling patterns, using the 
"master formula" in EqsJ3l|J36] combined with the simple (but accurate) parameterisa- 
tions of mixing angles discussed above, namely TBC in Eq|6]or TBC2 in Eq|4} We first 
recall the meaning of the ratios \zi\ in Eq|o"Tj namely that they correspond to the ratios 
of magnitudes of Dirac mass matrix elements, or equivalently, Yukawa couplings, and 
we repeat this result for convenience, 



\ z i\ 



l/l" 



Z 2 



\b\ 



*3 



(73) 



These ratios will be calculated in this subsection using Eqs,48 50 



Recall that the ratios of Yukawa couplings in the dominant column, z\ in Eq,48 



is 



independent of neutrino masses (and phases) and it only depends on the mixing angles 
and oscillation phase. In FigEq we see that \z\\ = 1 or \z\\ = 1/2 are possible for 
particular choices of physical phase 5 which depends on the precise values of the mixing 
angles. From the left panel of Figj2]we see that for TBC mixing we can achieve 



\Z\\ 



l/l 



(74) 
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Figure 2: The curves show the ratio of Yukawa couplings in the dominant column \z\\ = |e|/|/| as 
functions of the physical phase S evaluated using Eq|48j The left panel is for TBC mixing in Eq|6] while 
the right panel is for the TBC2 ansatz in Eq(4j both evaluated to second order in A. From the left 
panel we see that |e| = |/| is possible for 6 ~ ±7r/2 for TBC mixing. From the right panel we see that 
|e|/|/| « 1/2 is possible for S w the TBC2 ansatz in Eqji] 



for a particular choice of oscillation phase, S 
see that for TBC2 mixing we can achieve 



±7r/2. From the right panel of Figj2]we 



z i 



\e\ 

I7f 



i 

2 



(75) 



for a particular choice of oscillation phase close to zero, S w 0. 

These are quite interesting results, namely that simple ratios such as |e| = |/| or 
|e|/|/| = 1/2 are associated with simple values of oscillation phase, so we shall focus on 
these two examples. 

In the left panel of Figj3]we show the absolute magnitudes of the ratios \z2\ = \b\/\ a \ 
and \z 3 \ = \c\/\a\ as a function of the Majorana phase /3, for the case of TBC mixing, 
keeping the oscillation phase fixed at 5 — ir/2 corresponding to |e| = |/| in Eq|74j These 
ratios concern the second subdominant column of the Dirac mass matrix. As well as 
simple values of these ratios we are also interested in simple values of the physical phase 
differences 772 = 0b — <p e and 773 = C — 0/ and these are shown in the right panel of 
Figj3^is a function of the Majorana phase /3. 

By examining the results in Fig|3] we find some examples of simple ratios of Yukawa 
couplings in the subdominant column, which are shown together with the phases in 
Table HI Analogous results for 5 = —n/2 are also shown in Table [IJ where these results 
correspond to a symmetry under a change of the sign of all the phases 5, /3, i]2, ^3. This 
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Figure 3: Results for TBC mixing in EqJ6^ with S = n/2 and \zi\ — |e|/|/| = 1 as functions of the 
physical Majorana phase j3. The left panel shows the ratios of Dirac mass matrix elements j^l == H/I a l 
(solid) and \zs\ = \c\/\a\ (dashes). The right panel shows the physical phase differences ry 2 = 4>b — <j>e 
(solid) and 773 = (p c — (pf (dashes). Related results may be obtained for 5 = —ir/2 and \zi\ = |e|/|/| = 1. 
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Table 1: Examples of simple ratios of Yukawa couplings \z^\ = \b\/\a\, \z%\ = \c\/\a\ and corresponding 
approximate phases as estimated from Fig(3j The exact numerical PMNS parameters for selected cases 
will be considered later. 
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Figure 4: Results for TBC2 mixing in Eqpj] with 6 = and \z\_\ = |e|/|/| = 1/2 as functions of the 
physical Majorana phase /?. The left panel shows the ratios of Dirac mass matrix elements \z 2 \ — \b\/\a\ 
(solid) and \z^\ = |c|/|a| (dashes). The right panel shows the physical phase differences 772 = 4>b — 4>e 
(solid) and 773 = <j> e — <j>f (dashes). 



Mixing 


5 


w\ 


P 


N 


N 


m 


m 


TBC2 





0.5 





2.5 


1 








TBC2 





0.5 


±tt/4 


2 


1 


=F7r/6 





TBC2 





0.5 


±tt/2 


2 


2 





Tvr/2 


TBC2 





0.5 


±3tt/4 


2 


4.5 


Tvr/5 


T2tt/5 


TBC2 





0.5 


±7T 


1.5 


5.5 


Tvr/2 


Ttt/2 



Table 2: Examples of simple ratios of Yukawa couplings \z 2 \ = \b\/\a\, \z 3 \ — \c\/\a\ and corresponding 
approximate phases as estimated from Fig|3J The exact numerical PMNS parameters for selected cases 
will be considered later. 
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symmetry of the results can be understood by taking the complex conjugate of both 
sides of the "master formula" in Eqj29j We have verified this symmetry numerically. 

Similarly, in the left panel of Fig|4] we show the absolute magnitudes of the ratios 
I -22 1 — H/M an d | ^3 1 = |c|/|o| as a function of the Majorana phase 0, for the case of 
TBC2 mixing, keeping the oscillation phase fixed at 5 = corresponding to |e|/|/| = 1/2 
in Eqj75j Simple values of the physical phase differences rj2 = 4>b — 4>e and r] 3 = 4> c —(pf are 
shown in the right panel of Figj4]as a function of the Majorana phase /?. By examining 
the results in Figj4] we find some examples of simple ratios of Yukawa couplings in the 
subdominant column, which are shown together with the phases in Table [2] 

5 Brief Review of Indirect Models 

The results in the previous section can be used as the starting point to explore new types 
of indirect models. Indirect models are fully reviewed in [20J so here we only provide a 
brief review for completeness. 

Indirect models start from simple forms of the Yukawa columns which can be ob- 
tained from vacuum alignment due to some discrete non-Abelian family symmetry which 
is however completely broken. In the class of indirect models, the Klein symmetry of 
the neutrino mass matrix is not identified as a subgroup of the non-Abelian discrete 
family symmetry G. Models of this class are typically based on the type I seesaw mech- 
anism together with sequential dominance. Here, the main role of the family symmetry 



consists in relating the Yukawa couplings e, / as well as a, b, c of Eq. [47] by introducing 
flavon fields which transform as triplets under G (or anti-triplets if the representation 
is complex) and acquire special vacuum configurations. The directions of the flavon 
alignments are determined by the G symmetric operators of the flavon potential |20j . 

Continuing to work in the basis where both the charged leptons as well as the right- 
handed neutrinos are diagonal, the leptonic flavour structure is encoded in the Dirac 
neutrino Yukawa operator. The triplet (or anti-triplet) flavons in the neutrino sector 
denoted ^ atm , so i enter linearly as, 

C ~ ^LN atm H + ^fLN sol H + M atm A^ atm A atm + M sol A^ sol A sol + H.c (76) 

where A is a cut-off scale and H is the Standard Model Higgs doublet which is a singlet 
of G, while L are the three lepton doublets which form a triplet representation of G. 
The right-handed neutrinos A atm , A so i and the Higgs doublet H are assumed to be in 
the singlet representation of G. However to obtain only the above terms will require 
additional "shaping symmetries" which we do not specify here ( for examples which also 
include the charged lepton sector see the review [2D]). We obtain the Dirac neutrino 



Yukawa matrix by inserting the flavon VEVs into Eq. 76 Suppressing the dimensionless 
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couplings of the Dirac neutrino terms for notational clarity, we get 




,D ™D 



ni 



e.sol 



Y"v = m D = I e b I = j m? atm mf sol 1-^1 (0 a tm} 2 (0soi) 2 I • (77) 



/i,atm fi, sol 

m r,atm m r,sol 




The columns of the Dirac neutrino Yukawa matrix are therefore proportional to the 
vacuum alignments of the flavons fields a tm, 0soi- The effective Majorana operators 



of the light neutrinos can be derived from this using the seesaw formula of Eq. (16), 
yielding 

rv T ( (0atm) 1 (^atm) 7 (0sol) 1 (0sol) T ^ T c „ „ ,-,. 

^ ~ H~'^'^^ + ~'^r^^r ™' (78) 

which has the structure of EqjTTJ after the Higgs vev v is inserted. Note that the flavons 
enter the effective neutrino mass terms quadratically. 

In the class of indirect models, the PMNS mixing pattern thus becomes a question of 
the alignment vectors (0 a tm) and (0 so i). For instance, a neutrino mass matrix that gives 
rise to tri-bimaximal mixing can be obtained using the flavon alignments of constrained 
sequential dominance (CSD) 



(0 atm ) i r mi _ A.*. (<u 





A A Y MU V ; ' A 



Note that the resulting columns of the Dirac neutrino Yukawa matrix in Eq|77] are 
proportional to the columns of the unitary (in the present case TB ) mixing matrix with 
the physical phase difference (fib — (fie being related to the Majorana phase (3. Such a 
property of the Dirac neutrino Yukawa matrix is generally called form dominance [26J r\ 
Furthermore these alignments are left invariant under the action of the Klein symmetry 
generators, up to an irrelevant sign which drops out due to the quadratic appearance of 
each flavon. Since the family symmetry G does not contain the neutrino Klein symmetry, 
its primary role is then to explain the origin of these or similarly simple flavon alignments. 

It is clear that CSD is not realistic, primarily because it predicts TB mixing which 
involves a zero reactor angle. However the above example of CSD does illustrate the 
general strategy of obtaining simple ratios of Yukawa couplings within a particular 
column, including the texture zero d — 0, from vacuum alignment within the framework 
of indirect models using the techniques reviewed in [20] • 

An alternative approach to switching on the reactor angle without any "texture 
zeros" is to consider vacuum alignments which preserve the structure of CSD except that 
they allow d ^ 0, where such models are referred to as partially constrained sequential 



2 Exact form dominance implies vanishing leptogenesis |27j . This provides an independent motivation 
for deviating from CSD. 
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dominance (PCSD) [28]. The required vacuum alignment has been discussed within an 
explicit A 4 model in [29]. With the choice, 



i i4, \ (80) 

one obtains tri-bimaximal-reactor (TBR) mixing, which corresponds to TB mixing but 
with a non-zero reactor deviation parameter r and an oscillation phase 5 ~ 5q [28j 129] . 
For the choice r = X one arrives at TBC mixing [19]. However the main challenge for 
this approach is to explain why r = A, which is not easy to do (for an attempt see |19j). 
Another alternative approach is to consider particular alignments with two "texture 
zeros" . An example of this kind is, 

&f^ oc ( 1 J e * ^ oc [ 2 J e**> or ( J e** (81) 

These two alignment possibilities which involve two "texture zeros" were called CSD2 
[30] . Note that the above alignments involve the components of each of the flavon 
alignments (0 a tm) and (0 S oi) being relatively real. However the overall phase difference 
between the flavon alignments (pi, — <p e is physically significant and particular values of 
this phase difference corresponding to multiples of 7r/4 are preferred by the numerical 
fit resulting in maximal leptonic CP violation [30] • Unfortunately, as we have seen, it 
is not possible to obtain a large enough reactor angle for the case of a normal neutrino 
mass hierarchy with two texture zeros. Indeed the reactor angle predicted from CSD2 
was ^13 ~ 3 ^r leading to #13 ~ 5° — 6° which is too small [30J . The situation may be 
improved by considering charged lepton corrections [31] . but here we are considering the 
simple case where such corrections are completely negligible and in this case we cannot 
have two texture zeros in the two right-handed neutrino case of a normal neutrino mass 
hierarchy. 

Our approach here is to maintain a single "texture zero" d = which leads to 
a qualitative understanding of the reactor angle, 6*13 < m 2 /m 3 . However in order to 
provide a quantitative understanding of the lepton mixing angles, in the next section 
we shall consider particular types of vacuum alignment with (0 a tm)i = 0. 

6 Searching for New Types of Indirect Models 

In the previous subsection we briefly discussed the approaches which already exist in 
the literature, namely CSD, CSD2 and PCSD, which form the starting point of indirect 
models. In this section we shall use the results in this paper to explore new approaches 
which could form the starting point for new indirect models. The new models will 
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involve one single dominant texture zero, as in the case of CSD, and are unlike CSD2 
or PCSD which either involve two texture zeros or no texture zeros. However, whereas 
CSD predicts a zero reactor angle, the new approaches will predict a non-zero reactor 
angle consistent with experiment. 

To obtain realistic solutions which can accommodate the non-zero reactor angle, 
while maintaining the dominant texture zero d = (0 a tm)i = 0, we shall turn to the 
results of the previous section based on TBC and TBC2 mixing which are representative 
of the results from the latest global fits of leptonic mixing angles for the case of a normal 
neutrino mass hierarchy. From the point of view of Indirect Models we are looking for 
new examples where the magnitudes of the Yukawa couplings in a particular column 
occur in simple ratios. 

We emphasise that we only use TBC and TBC2 mixing to guide the search for simple 
alignments since they correspond to two phenomenologically acceptable mixing patterns. 
We thus allow simple alignments which only approximately satisfy TBC and TBC2 
mixing, and may be taken as the starting point for new models, leading to predictions 
for lepton mixing which do not fit exactly into either of these mixing patterns, but which 
we expect to be phenomenologically acceptable. 

We shall work in a phase basis where the dominant column involves Yukawa couplings 
which are real and positive, so all the phases occur in the second subdominant column. 
With this phase convention we can write, the Dirac mass (or Yukawa) matrix, previously 
given in Eq.68 together with Eq|77[ as, 





(82) 



where the phase (f) a is not physical and can always be adjusted to take any value by 
charged lepton re-phasing, while the two physical phases which are left invariant under 



charged lepton re-phasing, previously given in Eq.69 in this convention become 



V2 = 4>b, V3 = 4>c- (83) 

In this convention the flavon alignment (0 a tm) is purely real with a zero value in its first 
component and no phase difference between its second and third components, while the 
flavon alignment vector (<p B0 i) may in principle involve three phases, one for each compo- 
nent. Vacuum alignment techniques may lead to an arbitrary phase difference between 
different aligned components of the same flavon, as in the case of PCSD discussed in 



Eq.80 where the phase So is undetermined. However it is interesting to look for examples 
where 772 and 773 are are simply related, differing only by a simple rational multiple of n, 
since such simple relationships may be later understood in the context of some future 
theory. One particularly simple example would be where the components (</> so i)i, (0soi)2> 
(0soi)3 are relatively real, i.e. have the same phases up to ±7r (remembering that the 
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phase of (0 S oi)i is a free choice). An example of this kind of vacuum alignment is CSD2 



in Eq.81 A particularly simple example of this kind is, 



V = V2 = Vs- (84) 

In terms of general flavon models, the overall phases of the different flavon vevs (</> a tm) 
and (0soi) (i-e. the phases which are factored out of the alignments, for example the 



remaining phases which appear in CSD in Eq,79 or in CSD2 in Eq.81 ) are not determined 



since these overall phases depend on both arbitrary complex Yukawa coupling constants 
and complex flavon vevs. This is unfortunate since physical mixing angles and CP 
violating phases depend on this difference of overall phases, as in the example of CSD2 
discussed above. However, within special classes of models based on spontaneous CP 
violation, the Yukawa couplings are constrained to be real and in addition the phases 
of the overall flavon vevs may be determined using the vacuum alignment techniques 
discussed in [32] • Within such classes of models the overall phases of (0 a tm) and (0 so i) 
may be determined related in terms of simple rational multiples of n. Therefore we shall 
also be interested in such cases. 

From Table [T] we find the following phenomenological possibility, corresponding to 
TBC mixing with 5 = tc/2 and /3 = — 7r/4, where this case corresponds to rj = rj 2 = 




r/ 3 = — 27r/5, \z-\_\ = 1, \z 2 \ = 4, |z 3 | = 2. Using Eqs 82 and 83 we can see that this 
phenomenological example may be used as the basis for constructing a flavon model 
where the relative phases of the flavons are real, and the overall phases of the flavons 
take simple values, namely zero and — 2m/5, 

i^t oc | i ] = Aj ^ oc I 4 ] e - 2 -/ 5 = B, (85) 

where we have used the phase freedom in a to set this phase equal to —2tt/5. 

The alignments defined as A and B above are not sufficient to determine the effective 
neutrino mass matrix and hence the PMNS mixing parameters. Their normalisation is 
also required. In the case of form dominance, the alignments would be sufficient to 
determine the PMNS mixing parameters but form dominance is not satisfied. The 
effective neutrino mass matrix can be written as, 

m u = m a A T A + m b B T B = m a A T A + m a t p B T B (86) 

where m a and m& are real mass parameters which determine the physical neutrino masses 
m 3 and m 2 . For fixed alignments A, B, the ratio e u = m^jm a will also affect the PMNS 
mixing parameters. For the above alignment this ratio is given by, 




^=\-MA , (87) 
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where this ratio is related to physical parameters using Eqs j37]|40j For the case of TBC 
mixing with 5 = n/2 and /3 = — vr/4, with the mass ratio as in Eqj6lJ we find numerically 
that 



\b\ 



0.9159. 



This completely fixes the neutrino mass matrix up to an overall scale m a which does 
not affect the PMNS mixing parameters. Using the Mixing Parameter Tools (MPT) 
package associated with [33 J r\, we can check the predictions of this vacuum alignment, 
which may be compared to the original mixing pattern that was assumed to derive it 
(in brackets). For example, using the above alignment with e v = 0.058 we find using 
MPT, 7712/1713 = 0.17 together with, 
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34.2°, 0i3 = 9.2°, 6 23 = 40.9°, 5 = 107°, (5 



-83°, 



(89) 



which is in good agreement with the global fits. The discrepancy between Eqs|89] and 
the TBC mixing angles in Eq|5] arises due to the approximation in setting the ratios 
\zi\ and the phases r\i equal to simple values. We have checked that using the accurate 



values of the alignments, obtained from the master formula in Eq.29, with the ratios of 



Yukawa couplings calculated using Eqs|48j[49j[62]and the phases calculated using Eqj69 
together with Eqs, 70|7T][72" namely: 



A = 



with 



B = 



3.725e- i7r0 - 373 
2.028e-^°- 479 



(90) 



3.725 lei 



0.06045, 



(91) 



when used in EqJM] reproduce the TBC angles in Eq|5]to excellent accuracy. This gives 
us confidence that the master formula and the procedure that we followed to obtain the 
Yukawa ratios and the phases is correct. It also supports our approach of using the values 
obtained in Table [T] and Table [2] as the starting point for new simple but approximate 
alignments, which lead to phenomenologically viable possibilities which differ from either 
TBC or TBC2 mixing which were used to discover the new alignments, but which can 



subsequently be discarded. Thus the alignment in Eq.90 which reproduces TBC mixing 



is discarded in favour of the approximate but simpler alignment in Eq.85 which gives 



the different but phenomenologically viable mixing in Eq.89 



Note that the convention in 33 for the neutrino mass matrix m" differ from ours in Eq f 7 by an 



overall complex conjugation. Also the convention in 33J for the Majorana phases differs from ours in 



Eq 23 by a further complex conjugation. 
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As a second example from Table [T] we find the following phenomenological possibility, 
corresponding to approximate TBC mixing with 5 = tt/2 and /3 = 3tt/8, where this case 
corresponds to r\ = 7] 2 = 7/3 = tt/3, \z±\ = 1, \z 2 \ = 1, l^l = 3. This suggests the 
alignments, in our phase convention, 

A = j 1 ] , B = I 1 I e in/3 . (92) 

In a realistic model, the alignments may be predicted by some family symmetry for 
example, but not the value of e u which will depend on right-handed neutrino masses 
and unknown Yukawa couplings. Therefore we shall regard e u as a free parameter which 
may be varied to give different the correct ratio of physical neutrino masses 7712/1713. For 
example, using the above alignment with t v = 0.1 we find using MPT, 7712/771-3 = 0.17 
together with, 

# 12 = 34.3°, 13 = 8.6°, # 2 3 = 44.3 , 5 = 93°, f3 = 72°, (93) 

which is in good agreement with the global fits. 

As a third example from Table [T] we find the following phenomenological possibility, 
corresponding to approximate TBC mixing with 5 = it/2 and /3 = —3tt/8, where this 
case corresponds to r\ = r/ 2 = r/ 3 = — tt/3, \z±\ = 1, \z 2 \ = 3, \z 3 \ = 1. This suggests the 
alignments, in our phase convention, 

A = j 1 ) , B = I 3 I e~ OT/3 . (94) 

For example, using the above alignment with e u = 0.1 we find using MPT, 7712/7713 = 0.17 
together with, 

fl 12 = 34.3°, #13 = 8.6°, # 2 3 = 45.7°, 5 = 87°, (3 = -72°, (95) 

which again is in good agreement with the global fits. 

We emphasise that, with the alignments including the phase 7] fixed, the neutrino 



mass matrix is completely determined by only two parameters in Eq|86| namely an 
overall mass scale m a , which may be taken to fix the atmospheric neutrino mass 777.3 = 
0.048 — 0.051 eV, the ratio of input masses e„, which may be taken to fix the solar to 
atmospheric neutrino mass ratio 7772/7773 = 0.17 — 0.18. In particular the entire PMNS 
mixing matrix and all the parameters therein are then predicted as a function of 7772/7773 
controlled by the only remaining parameter e„. In Table [3] we show the predictions for 



the first example above in Eq 85 as a function of e„ and hence 7712/7773. In Table 0we 
show the predictions for the second example above in Eq|92]as a function of e u and hence 
7712/777,3. In Table [5] we show the predictions for the third example above in EqJ94]as a 
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e v 


m 2 /m 3 


#12 


#13 


#23 


5 


P 


0.057 


0.166 


34.2° 


9.0° 


40.8° 


107° 


-84° 


0.058 


0.170 


34.2° 


9.2° 


40.9° 


107° 


-83° 


0.059 


0.174 


34.1° 


9.4° 


41.0° 


106° 


-82° 


0.060 


0.177 


34.1° 


9.6° 


41.1° 


105° 


-80° 


0.061 


0.181 


34.1° 


9.7° 


41.3° 


104° 


-79° 



Table 3: The predictions for PMNS parameters and 7712/7123 arising from Eq 85 as a function of e„. 
Note that these predictions assume r\ = —2ir/5. Identical results are obtained for r\ = 2ir/5 with the 
phases j3 and 8 changed in sign. 



e v 


m 2 /m 3 


#12 


#13 


#23 


5 


P 


0.098 


0.168 


34.4° 


8.4° 


44.4° 


92° 


-73° 


0.100 


0.171 


34.3° 


8.6° 


44.3° 


93° 


-72° 


0.102 


0.173 


34.3° 


8.75° 


44.1° 


94° 


-71° 


0.104 


0.177 


34.3° 


8.9° 


44.0° 


94° 


-70° 


0.106 


0.179 


34.2° 


9.1° 


43.8° 


95° 


-69° 



Table 4: The predictions for PMNS parameters and 772,2/7713 arising from CSD3 in Eq 92 as a function 
of e„. Note that these predictions assume 77 = it/ 3. Identical results are obtained for 77 = — n/3 
with the phases j3 and 6 changed in sign. These angle predictions are approximately equal to the 
Tri-bimaximal-Cabibbo mixing angle values #12 = 35.26°, #13 = 9.15°, 623 = 45° to within one degree. 



function of e v and hence 1712/7113. The results in Table fi] and Table p^ only differ in the 
atmospheric angle and the oscillation phase S, which are correlated via an atmospheric 
sum rule discussed later in Eq ]153[ with both sets of results approximated by TBC 
mixing to within an accuracy of one degree, but with the phases predicted. 

We remark that an accuracy of one degree in the angles is all that can be expected 
due to purely theoretical corrections in a realistic model due to renormalisation group 
running [34] and canonical normalisation corrections [35]. In addition, there may be 
small contributions from a heavy third right-handed neutrino [36] which can affect the 
results. 



7 Leading Order Analytic Results 



The above results are numerical, based on the exact analytic master formula in Eq.29 



In order to get some feeling for the results it is worth taking a look at the well known 
but approximate analytic expressions for the masses and mixing parameters in terms 
of the Dirac mass matrix elements and right-handed neutrino masses, bearing in mind 
that these results may have large corrections. The neutrino masses in terms of the 



25 



e v 


m 2 /m 3 


#12 


#13 


#23 


5 


P 


0.098 


0.168 


34.4° 


8.4° 


45.5° 


88° 


-73° 


0.100 


0.171 


34.3° 


8.6° 


45.7° 


87° 


-72° 


0.102 


0.173 


34.3° 


8.75° 


45.9° 


86° 


-71° 


0.104 


0.177 


34.3° 


8.9° 


46.0° 


86° 


-70° 


0.106 


0.179 


34.2° 


9.1° 


46.2° 


85° 


-69° 



Table 5: The predictions for PMNS parameters and 7712/7713 arising from CSD3 in Eq 94 
of e 



as a function 
it/3. Identical results are obtained for r\ = tt/3 



Note that these predictions assume 77 = 
with the phases j3 and 6 changed in sign. These angle predictions are approximately equal to the 
Tri-bimaximal-Cabibbo mixing angle values #12 = 35.26°, #13 = 9.15°, #23 = 45° to within one degree. 



parameters in Eqjlllare given to leading order in 1712/7113 by [7J HP] , 

mi = 

\a\ 2 
m * * X4 2 

\e\ 2 + \f\ 2 



777,3 



Y 



The neutrino mixing angles are given to leading order in 777.2/717.3 by [7J HP] , 



tan# 



23 



l/l 



tan# 



12 



9 



13 



c 23 |6| cos(> 6 ) - s 23 |c| cosOJ 



[|e| 2 + |/| 2 ]3/2X 



(96) 
(97) 

(98) 



(99) 
(100) 
(101) 



where we have written some (but not all) complex Yukawa couplings asi= \x\e l( ^ x . The 
Dirac CP violating oscillation phase 5 is fixed to give a real angle #12 by, 



c 23 \b\ sm(4> b ) « s 23 |c|sin( 



(102) 



where 

4>b = 4>b - 4>a - 4> + 8, 

4> c = 4> c - 4> a + 4> e - 4>f - 4> + 5 

The phase is fixed to give a real angle 6*13 by [7], 



(103) 



(104) 
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where 



arg(e*6 + f*c). 



Eq ]105| may be expressed as 

, , _ \b\s 23 s 2 + \c\c 23 s 3 

tcLll (j) ~ —— : — : . 

|0|S 2 3C 2 + \C\C23C3 

Inserting <fi in Eq, 104 into Eqs, 102|l03 , we obtain, 

tan(0 + ^ 1^3*' -H***3 
-\b\C23C2 + \c\S23C3 

where we have written Sj = sinr/j, q = cos rji where 



(105) 



(106) 



(107) 



(108) 



are the phases introduced in Eq.69 which are invariant under a charged lepton phase 
transformation. The reason that the see-saw parameters only involve two invariant 
phases 772 , V3 rather than the usual six is due to the two-right handed neutrino assump- 
tion, which removes three phases, together with the assumption of a dominant texture 
zero, which removes another phase. 

From the above approximate leading order results we find the following alternative 
useful expressions for the mixing angles, 



tan# 



23 



COt0 



12 



'13 



\Zi\ 
C 23 |2 2 |C0S I 1)2 



S 23 |2 3 |C0S 773 



P 



77_V 

m 3 " 



s? 2 C23|N + Ntanfl 2 3e J (")| 



(109) 
(110) 

(111) 



where t; 2 and 1)3 are given in Eq.69 and 



are given in Eq,73 The above approximate 

In fact 



leading order results should be compared to the exact results in Eqs|48]j50 
they may be derived by inverting 



48]|50] This is non-trivial to do exactly, but to 

using the standard PDG 
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leading order in m 2 /?7i3 it leads to the results in Eqs. 109 
expressions for the mixing elements. We emphasise that these leading order results 
are very crude, for example we have seen from FigfJJthat for close to maximal mixing 
Z\ w 1 has a 50% variation depending on the phase 5. Indeed maximal atmospheric 
mixing prefers 5 ~ n/2. This is due to next to leading order effects beyond the simple 
approximations above. However the leading order results do provide some understanding 
of the successful alignments discussed in the previous subsection as we now discuss. 
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The alignments discussed in the previous subsection have 77 = 7/2 = 773 and lead to 
approximate maximal mixing with \zi\ = 1 and tan #23 ~ 1. In this case the approximate 
formulas above become, 

cot9 12 « ^cosh-0(|^ 2 |-NI) (112) 

Ois - %^(M + M)- ( 113 ) 

These approximate results show that the solar angle depends on the difference between 
I Z2 1 and I213 1 and the phase difference 77 — | . Phenomenologically we must arrange to 
obtain cot 9 12 ~ v^ leading to the approximate condition for the difference, 

N-N « — ; «T (114) 

COS{1] - I J 

Since s^ 2 ~ 1/3 the reactor angle is approximately given by, 

^ " 371 S (N + N) ' (U5) 

Since phenomenologically 7122/1113 ps A/a/2 and 6*13 ~ A/\/2, we obtain the approximate 
condition for the sum, 

|z 2 | + |z 3 |w3>/2. (116) 

In the examples where r? = r? 2 = r? 3 , there is only a single see-saw phase 77 and so 
all physical phases must be related to it. In particular we find, using the above results 
and other approximate results for the phases in [7J [10] that, for |z 2 | > \zs\, </> = V and 
5 ~ —2r) and /3 ~ 2rj. Thus the Majorana phase (3 is simply related to the oscillation 
phase 5 according to 5 ~ — /3. From Eq ]114| we find in these cases, 

\z 2 \-\z 3 \ w 2. (117) 

Solving Eqsjll6| and |117| we find, 



|z 2 |«-p + lw3.1, |* 3 | « -p - 1 « 1.1- (118) 

1 V2 ' V2 

These results are well respected by the third example in EqJ94] corresponding to \z\\ — 1 



I Z2 1 = 3, I £3 1 = 1. For the second example analogous results to those in Eq. 118 may 



be derived but with |z 2 | and \z^\ being interchanged, which is in good agreement with 



\z%\ = 1, I ^3 1 = 3. The first example in Eq,85 with \z±\ = 1, |z 2 | — 4, \zs\ = 2 does 



not obey the estimate in Eq |118 very well. Nevertheless it leads to good agreement 



with experiment, which highlights the limitations of the leading order results, and the 
importance of the exact master formula. 
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8 Link between Leptogenesis and PMNS phases 

In the two right-handed neutrino model with dominant texture zero, that we are consid- 



ering in this paper, Eqj83] shows that there are only two see-saw phases that are relevant, 
namely 772 and 773. The low energy PMNS phases 5 and /3 are therefore (complicated) 
functions of these two see-saw phases, and the other see-saw parameters, with the exact 



relation obtained via the master formula in Eq. 29 The phases relevant for leptogenesis 
can also only depend on the two see-saw phases 772 and 773. This implies that the phases 
which control leptogenesis must be related to the two physical low energy PMNS phases 
S and j3. This means that by studying CP violation in the laboratory one is in fact also 
probing CP violation in the early universe which was relevant for matter-antimatter 
asymmetry. 

The Leptogenesis-PMNS link for the two right-handed neutrino model with dominant 
texture zero been discussed some time ago in the literature, first for flavour independent 
leptogenesis in [10], then later for flavour dependent leptogenesis in [37J. In the flavour 
independent leptogenesis case [10J the lepton asymmetry is proportional to sin 2<ft where 



the relevant leptogenesis phase is <p was defined in Eq. 105 In the flavour dependent 
leptogenesis case [37] muon-type lepton asymmetry is proportional to sin(0 + 772), with 
tau-type lepton asymmetry being proportional to sin(0 + 773) and electron-type lepton 
asymmetry being identically zero due to the texture zero. 

Estimates of the baryon asymmetry Yb, including the effects of washout, for the case 
where the lightest right-handed neutrino dominates the production of baryon asymmetry 
(so-called iVi leptogenesis) for the flavour-dependent case were made in [37] , in particular 
for the above case of the two right-handed neutrino model with a dominant texture zero. 
When the lightest right-handed neutrino is identified with iV atm , one finds [37], assuming 
|e| ~ |/| which leads to equal washouts r/ M ~ r\ T 



IT- 



Y B oc |b| sin(0 + r] 2 ) + \c\ sin(0 + 773). (119) 

When the lightest right-handed neutrino is identified with iV sol , one finds [37], assuming 
M ~ l/l) with unequal washouts in this case, 

Y B oc -r]^\b\ sin(0 + 772) - 77 r |c| sin(0 + 773). (120) 

The above estimates assume that the lightest right-handed neutrino dominates the pro- 
duction of baryon asymmetry (iVi leptogenesis). In certain situations the heavier right- 
handed neutrino can dominate (N 2 leptogenesis) as was discussed for the two right- 
handed neutrino model, especially for the case where the lightest right-handed neutrino 
is identified with iV atm [38]. In this case the situation is more involved, and it is not 
possible to give simple analytic estimates for the baryon asymmetry as above. How- 
ever we remark that for N 2 leptogenesis the flavour dependent lepton asymmetries are 
proportional to the same phase factors as for Ni leptogenesis but have opposite sign. 
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For the examples where rj = r]2 = 773 the baryon asymmetry Yb is clearly proportional 



to a common factor sin(0 + 77). From Eq. 105 we see that in this case 4> = rj which is an 



exact result. Hence the baryon asymmetry is proportional to sin(2r/) or sin(2</>) which 
makes sense since 77 is the only remaining see-saw phase in this case. The physical low 
energy CP violating phases and 5 must vanish in the limit that 77 — > since 77 is the only 
high energy source of CP violation at the see-saw scale. The exact relationship between 
(3, 5 and 77 provided via the master formula is rather complicated, however the leading 
order estimates for the phases derived from pj [TU] are 5 ~ — 2?7 and (3 rj 2?7, which are 



in fair agreement with the numerical results in Eq.95 Thus, approximately, we have in 
this case that the lepton asymmetry is proportional to sin S, which is maximal in this 
case. The main message is that, by studying CP violation in oscillation experiments, 
one is in fact also probing the source of CP violation in the early universe which was 
relevant for matter-antimatter asymmetry. 

9 A Minimal Predictive See-Saw Model 

9.1 CSD3 and vacuum alignment 

In discussing the implications of our results for model building, we have been guided 
by the idea that the most attractive flavon alignments are the simplest ones, since they 
may eventually arise from some realistic model. Although we have tabulated several 
relatively simple possibilities, we have only explored three examples from Table [T] in 
detail. These three examples all satisfy the conditions that 77 = 772 = 773 so that the 
components of a given flavon all carry the same phase and lead to approximate maximal 
mixing with \z±\ = 1 which means that the atmospheric flavon alignment is particularly 



simple and retains the CSD form in Eq. 79 However the price of the non-zero reactor 



angle (at least in our minimal approach where we maintain the zero in the atmospheric 
vacuum alignment, and do not invoke charged lepton corrections) is that the solar flavon 



alignment must differ from the CSD alignment shown in Eqf79| and in general will be 
more complicated. The question then arises as to whether such alignments may be 
achieved theoretically from vacuum alignment methods. 

In this subsection we shall discuss two of the successful flavon alignments which 
we have discovered in our general analysis, and show how they may be derived from 
vacuum alignment methods. The two examples we shall consider (which correspond to 
the second and third examples discussed earlier in this section) are summarised below, 

'*•*> oc [ 1 I e'*', %^(x|l|e*or I 3 1 1121; 



A \i/ A W V, 

These two new alignment possibilities shall both be referred to as CSD3, since they 
both represent new possibilities not considered before in the literature. As in the case of 
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CSD2 [30] the above alignments involve the components of each of the flavon alignments 
(0atm) and (0soi) being relatively real within each flavon, although the overall phase of 
that flavon will in general be non-zero. However the overall phase difference between the 
flavon alignments r\ = <pb — 4> e is physically significant and we have seen that particular 
values of this phase difference are preferred. In a particular phase convention, the 



preferred values are displayed in Eqs.92 and 94 

The vacuum alignments associated with TB mixing have been very well studied. 
Here we shall focus on the family symmetry A4 as it is the smallest non-Abelian finite 
group with an irreducible triplet representation. We apply the A± basis in which the 
triplets are explicitly real as given for example in [39]. We shall work in the basis of 
[29] where, denoting a general A 4 triplet as c = (c\, c 2 , cz) T and defining to = e 27 ™/ 3 , the 
product rules can be summarised as 



c' 



r=0 




Y,(cA + u~ r c 2 c' 2 + u r c 3 c' 3 ) + c 3 c[ + Cl 4 , (122) 



corresponding to two triplets and the sum of the three one-dimensional irreducible rep- 
resentations l r , with 1 = 1 being the trivial singlet. 

Following the methods of [29] it is straightforward to obtain the vacuum alignments 
for charged lepton flavon alignments suitable for a diagonal charged lepton mass matrix. 
The charged lepton flavon alignments used to generate a diagonal charged lepton mass 
matrix are obtained from the renormalisable superpotential [29J, 



rrrflavon,! 
A4 



r^ 



A e ip e ip e + AyWiMp + A T ip T ip T + O efl tp e tpv + O eT (p e (p T + O^if^r .(123) 

The triplet driving fields A e ^^ T give rise to flavon alignments ((p e ,n,r) with two zero 
components, and the singlet driving fields Oij require orthogonality among the three 
flavon VEVs so that we arrive at the vacuum structure 129 1, 



(<p e ) = v e \0] , M = vjl\ , ((p T ) = v T [0] . (124) 




Of more interest to us in this paper are the new neutrino flavon alignments. The starting 
point for the discussion is the usual standard TB neutrino flavon alignments proportional 
to the respective columns of the TB mixing matrix, 



(W> = v Vl -1 , ((fu 2 ) = v U2 \ 1 , (^ 3 ) = v U3 1 . (125) 
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In the remainder of this subsection we shall show how to obtain the neutrino flavon 
alignments including the new alignment, 

V 



(<Pu t ) 



"4 



(126) 



which is an example of a CSD3 solar flavon alignment in Eq ]121[ In constructing a 



CSD3 model in the notation of Eq|121 we shall identify a t m 



<r 



Ki ■ 



and 



^sol 



<r 



VA- 



The renormalisable superpotential involving the driving fields necessary for aligning the 
neutrino-type flavons is given as^ 



W 



flavon,!/ 

A 4 



Au 2 (gi<Pu 2 <Pv 2 + g2<Pva€va) (127) 

+ O eV3 g 3 (p e (p V3 + U2U . a g^u 2 fu 3 + UlU2 g 5 ip Ul ip U2 + UlU3 g 6 ip Ul ip U3 
+ O fiU5 g 7 (p IM (p V6 + U2Ub gs(p V2 ip Vb + O^ggippipvg + U5VQ g 1Q ^ U5 (p v& 

+ V6Vi guLp vr> ip Vi + UlVi gi2^ Ul ^u^ 

where A V2 is a triplet driving field and O^ are singlet driving fields whose F-terms 
lead to orthogonality relations between the accompanying flavon fields. Here gi are 



dimensionless coupling constants. The first line of Eq. (127) produces the vacuum 
alignment (<p U2 ) oc (1, 1, — 1) T of Eq. ( |125 ) as can be seen from the F-term conditions^ 

2<?i (^) 3 (^ 2 )i + <&<&*> (^)2 = . (128) 



The first two terms in the second line of Eq. ( 127 ) give rise to orthogonality conditions 



which uniquely fix the alignment (fi^) oc (0, 1, 1) of Eq. (125) 



(Vef ■ (<fu 3 ) = (<fu 2 ) T ■ ((fu 3 ) 







-> 



((fu 3 ) OC 




(129) 



The last two terms in the second line of Eq. (|127|) give rise to orthogonality conditions 
which uniquely fix the alignment (<f vi ) oc (2, 



(ip Ul ) T ■ ((p U2 ) = ( Vui ) T ■ (ip U3 ) 








(130) 



4 The alternative CSD3 alignment {<p Ui ) °c (1, 1, 3) T may readily be obtained by replacing ip^ — > ip r 



127 



in Eq 

5 We remark that the general alignment derived from these F-term conditions is (</V 2 ) c>c 
(±1,±1,±1) T . One can, however, show that all of them are equivalent up to phase redefinitions. 
Note that (1,1,-1) is related to permutations of the minus sign as well as to (—1,-1,-1) by A4 trans- 
formations. The other four choices can be obtained from these by simply multiplying an overall phase 
(which would also change the sign of the £„ 2 vev.) 
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The terms in the third line of Eq. (127) give rise to orthogonality conditions which 
uniquely fix the alignments of the auxiliary flavon fields (p V5 and (p ue , 



(¥V) T • (<Pvs) = (^ 2 ) T • (Vl*) = ~* fans) « I 



(131) 



(^) T ■ (^ 6 ) = ( < Pu 5 ) T • (<^ 6 > =° -> fam) OC j 



(132) 



The neutrino-type flavon of interest labelled as y? V4 gets aligned by the remaining 
terms in the fourth line of Eq. (127). A vanishing F-term of the driving field U6Ui 
requires 




(W> = I "2 I , (133) 

where n\ and n 2 are independent parameters. Finally a vanishing F-term of the driving 



field UlUA leads to the desired alignment in Eq |126 







— >■ 




(134) 



So far we have only shown how to align the flavon vevs and have not enforced them 
to be non-zero. In order to do this we shall introduce the additional non-renormalisable 
superpotential terms which include, 



AW 



flavon/ 

A 4 



E j((^-Vi)Pi-M 3 ) + j(pf-M 3 ), 



l=e,fi,r 



aw m ~ 



Ei(W'^.- m3 ) + t«- m3 



AW 



flavon, £„ 

A 4 



t=l 
P 



{it - m a 



A v^ 2 



(135) 

(136) 
(137) 



where, as in [32], the driving singlet fields P denote linear combinations of identical 
singlets and we have introduced explicit masses M to drive the non-zero vevs, as well as 
the messenger scales denoted as A. For example the A± singlet £„ 2 has a vev of order M 



driven by the F-term of the singlet P in Eq, 137 We have also introduced A^ singlets 
Pi and p Ui whose vevs are driven by the F-terms of the singlets P in the second terms 



in Eqs, 135 and 136 These singlet vevs enter the first terms in Eqs. 135 and 136 which 
drive the vevs of the triplet flavons. 
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The flavons and driving fields introduced in this subsection transform under Z\ x Z^ 



symmetries whose purpose is to allow only the terms in Eqs ]123[ |127 and 135 137 and 



forbid all other terms. Any superfield with a single subscript V{ transforms under Z^ 
as u 2 and is a singlet under all other discrete symmetries. Any superfield with a single 
subscript I transforms under Z\ as u and is a singlet under all other discrete symmetries. 
The orthogonality driving superfields Oij with two subscripts transform under Z\ x Z^ 



in such a way as to allow the terms in Eqs. 123 127 For example the 0\ Vi driving fields 
transform under Z\ x ZQ as (u> 2 , u). In addition driving superfields are assigned a charge 
of two while flavon superfields have zero charge under a U(1)r symmetry. 

9.2 A predictive A4 x Z\ Q model of leptons 

In this subsection we outline a supersymmetric (SUSY) v4 4 model of leptons with CSD3 
along the lines of the A 4 models of leptons discussed in [301 EH]- The basic idea is that the 
three families of lepton doublets L form a triplet of A4 while the right-handed charged 
leptons e c ,/i c ,r c , right-handed neutrinos Na_t m ,N so i and the two Higgs doublets Hi,H 2 
required by SUSY are all singlets of A4. In addition the model employs an additional 
discrete symmetry Z\ x Z^- introduced in the previous subsection together with a Z^ 
family symmetry in order to account for the charged lepton mass hierarchy. The full 
discrete group of the model is A4 x Z\° which has order 42. 

In Table [6] we have displayed the symmetries and superfields relevant for the Yukawa 
sector and have written Z| tm = Z^ 3 and Zf l = Z3 4 . The transformation properties 
of the remaining superfields under Z\ x Z^ was discussed at the end of the previous 
subsection and are consistent with the charges shown in Table [6} 

The charged lepton sector of the model employs the A± triplet flavons (p e , (p^, <p T 
whose alignment was discussed in the previous subsection. With the lepton symme- 
tries in the upper left of Table [6] we may enforce the following charged lepton Yukawa 
superpotential at leading order 

W Yuk ~ iffifov ■ L)t c + ^MTifo, • L)// + Y^H^e ■ L)e c , (138) 

which give the charged lepton Yukawa couplings after the flavons develop their vevs. A is 
a generic messenger mass scale, but in a renormalisable model the messengers scales may 
differ. The charged lepton symmetries include three lepton flavour symmetries Z^' T 
under which ip e ,ip fJi ,if T and e c ,/i c ,r c transform respectively as uj and u 2 , together with 
a lepton family symmetry Z® under which e c ,/i c ,r c transform as u,u 2 ,l respectively 
(where u = e l27r ^ 3 ) with the family symmetry breaking flavon 9 transforming as to and 
otherwise being a singlet under all other symmetries. H\ and L and all other fields are 
singlets under Z^' r and Zf. With these charge assignments the higher order corrections 
are very suppressed. 
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Table 6: Lepton, Higgs and flavon superfields and how they transform under the symmetries relevant 
for the Yukawa sector of the model. The only non-trivial charged lepton charges are in the upper left 
of the Table and the only non-trivial neutrino charges in the lower right of the Table. Note that the 
only the lepton doublets L and A4 symmetry, are common to both charged lepton and neutrino sectors 
and are given near the central column and row. The Standard Model gauge symmetries and U(1)r 
symmetry, under which all the leptons have a charge of unity while the Higgs and flavons have zero 
charge, are not shown in the Table. 

The charged lepton Yukawa matrix is diagonal at leading order due to the alignment 



of the charged lepton-type flavons in Eq, 124 (where the driving fields responsible for the 



alignment in Eq. 123 absorb the charges under the newly introduced symmetries Z 



e,n,r 



and Z|) and has the form, 



Y e = diag(y e , j/ M , y T ) ~ diag(e 2 ,e, l)y T 



(139) 



where we choose e ~ (#)/A ~ A 2 in order to generate the correct order of magnitude 
charged lepton mass hierarchy, with precise charged lepton masses also dependent on 
order one coefficients which we have suppressed here. 

The neutrino sector of the model exploits the A 4 triplet flavons at m = <Pv 3 , and 
sol = (p U4 whose alignment was discussed in the previous subsection. We shall consider 
one of the types of CSD3 corresponding to, 
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(140) 



With the neutrino symmetries in the lower right part of Table [6] we may enforce the 
following leading order neutrino Yukawa superpotential 



W Yuk ~ X # 2 (0at m 



£)iV atm + ^# 2 (0 sol -L)iV sol . 



fi4r 
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Again the higher order corrections are completely negligible. Inserting the vacuum 
alignments in Eq jl40| the neutrino Yukawa matrix is of the form 



Y v = e 36 . (142) 




The parameters e and b can be determined from the parameters in the superpotential. 
As is typical in models of this kind [301 129] , the RH neutrinos have no mass terms 
at the renormalisable level, but they become massive after some A 4 singlet flavons £ atm 
and £ so i develop their vevs due to the renormalisable superpotential, 

W*~UX tm + £ sol iV s 2 ol . (143) 

When the right-handed neutrino flavons develop their vevs (£ a tm) ~ Ma together with 
(£soi) ~ Mb> then the RH neutrino mass matrix is diagonal as required, 

*■-(?£)■ (144) 

To ensure that the mixed terms are absent at renormalisable order we have imposed a 
right-handed neutrino flavour symmetry Z^ tm under which iV a t m and £ atm both transform 
as to (where to = e* 27r//3 ) while a tm transforms as co 2 , with all other fields being singlets. 
We have also imposed a similar symmetry Zf l under which the "solar" fields transform 
in an analogous way. We remark that these charge assignments are consistent with the 



flavon superpotential in Eq 127, where we identify atm = (p U3 , and so i = ip V4 , with 
suitable charges assigned to the driving fields. 

The above charge assignments allow higher order non-renormalisable mixed terms 
such as 

AW R ~ ^(0 a tm-0 S ol)^ atm AT sol , (145) 

which contribute off-diagonal terms to the right-handed neutrino mass matrix of a mag- 
nitude which depends on the absolute scale of the flavon vevs (</> a tm) and (0 so i) compared 
to (£ a tm) and (^soi)- If all flavon vevs and messenger scales are set equal then these terms 



are suppressed by e ~ A 2 according to the estimate below Eq, 139 however they may be 
even more suppressed. We shall ignore the contribution of such off-diagonal mass terms 
in the following. 

The effective neutrino mass matrix after implementing the see-saw mechanism can 
be written from Eqs J86|l40fT44 , 



'0 (T 

in--. m a A T A + m a e 2iv e„B T B = m a \0 1 1 ] + m a e 2ir >e v ( 3 9 3| (146) 

,0 11, 
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where m a and e„ are real mass parameters which determine the physical neutrino masses 
777,3 an d "72 and the phase difference 77 was defined in Eq,84, Thus the neutrino mass 
matrix is determined by three input parameters: m a , e u and 77. These three input 
parameters determine nine physical observables consisting of the PMNS mixing matrix 
plus the three neutrino masses, making this a highly predictive scheme. We saw earlier 
that a relative phase difference 77 = — 7r/3, which translates into a Dirac CP phase 
5 ph 7r/2 gives a successful set of mixing parameters. 

Such a phase difference 77 = — 7r/3 could be achieved in the context of spontaneous 
CP violation from discrete symmetries as discussed in [32], and we shall follow the 
strategy outlined there. The basic idea is to impose CP conservation on the theory so 
that all couplings and masses are real. Note that the A4 assignments in Table [6] do not 
involve the complex singlets 1', 1" or any complex Clebsch-Gordan coefficients so that 
the definition of CP is straightforward in this model and hence CP may be defined in 
different ways which are equivalent for our purposes (see [32] for a discussion of this 
point). The CP symmetry is broken in a discrete way by the form of the superpotential 
terms. We shall follow [32J and suppose that the flavon vevs (0 a t m ) and (0 so i) to be real 



with the phase 77 in Eq. 146 originating from the solar right-handed neutrino mass due 
to the flavon vev (£ so i) ~ M#e 2l7r//3 having a complex phase of 2n/3, while the flavon 
vev (£atm) ~ M^ is real and positive. This can be arranged if the right-handed neutrino 
flavon vevs arise from the superpotential, 
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(147) 



where, as in [32], the driving singlet fields P,P' denote linear combinations of identical 
singlets and all couplings and masses are real due to CP conservation. The F-term 
conditions from Eq jl47| are, 
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These are satisfied by («£ atm ) = ((AM 2 ) 1 / 3 ] and (£ sol ) = |(A , M ,2 ) 1 /3| e 2^/3 where we 
arbitrarily select the phases to be zero and 27r/3 from amongst a discrete set of possible 
choices in each case. More generally we require a phase difference of 27r/3 since the 
overall phase is not physically relevant, which would happen one in three times by 
chance. In the basis where the right-handed neutrino masses are real and positive this 
is equivalent to r] = — 7r/3 in Eq, 146 due to the see-saw mechanism. 



Similarly the flavons appearing in Eqs, 135 137 each have a discrete choice of phases 
(0, 27r/3,47r/3). The charged lepton flavons ipi may take any phases since such phases 
are unphysical. In fact the only physically significant flavon phases from the previous 
subsection are those of atm = (p U3 
to maintain 77 = —it/3 in Eq, 146 
chance. 



and so i = ip U4 whose phases are selected to be equal 
As before, this would occur one in three times by 
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In Table [5] we show the predictions for this model corresponding to Eqj94] fixing 
7] = —tt/3 as a function of e v and hence m 2 /m3. With rj = — n/3, the PMNS parameters 
are all accurately predicted and only depend on the physical neutrino mass ratio mijvn^ 
with mi = 0. This level of predictiveness for CSD3 is equal to that for CSD2. However 
in the case of CSD2 the reactor angle was wrongly predicted, whereas CSD3 reproduces 
the observed value very well. 

As in the case of CSD2, the neutrino mass matrix implies the TM4 mixing form |40] 
where the first column of the PMNS matrix is proportional to (2,-1, 1) T . The reason is 



simply that (<f Ul ) ^ (2, — 1, 1) T is an eigenvector of m u in Eq 146 with a zero eigenvalue 
corresponding to the first neutrino mass m\ being zero. The reason for this is that m v 
in Eq, 146 is a sum of two terms, the first being proportional to v4v4 T oc (ip v , i )((p Va ) T 



and the second being proportional to BB T oc ({p U4 )({p Ui ) T . Since (<p Ul ) oc (2, — 1, 1) T is 
orthogonal to both (<f V3 ) and (<f V4 ) it is then clearly annihilated by the neutrino mass 
matrix, i.e. it is an eigenvector with zero eigenvalue. Therefore we immediately expect 



m" in Eq. 146 to be diagonalised by the TMi mixing matrix [40j where the first column 
is proportional to (w) °^ (2, — 1, 1) T - Therefore we already know that CSD3 must lead 
to TMi mixing exactly to all orders according to this general argument. 

Exact TMi mixing angle and phase relations are obtained by equating moduli of 
PMNS elements to those of the first column of the TB mixing matrix: 

C12C13 = J-, (149) 

IC23S12 + si 3 S23Ci2e ?<5 | = —=, (150) 

V6 

IS23S12 - Sl3C 2 3Cl2e j5 | = —7=. (151) 



From Eq jl49| we see that TMi mixing approximately preserves the successful TB mixing 
for the solar mixing angle 9 12 ~ 35° as the correction due to a non-zero but relatively 
small reactor angle is of second order. While general TMi mixing involves an un- 
determined reactor angle #13, we emphasise that CSD3 fixes this reactor angle. The 
approximate leading order result from Eq |115| is 

which explains why the bound #13 < 777,2/^3 is approximately saturated. However the 
leading order results are not highly accurate and numerically the prediction for the 
reactor angle depends on the phase n. For rj = 7r/3 the reactor angle is in the correct 
range as shown in Table [5j 



In an approximate linear form, the relations in Eq. 149 151 imply the atmospheric 
sum rule relation, 

# 23 ^45° + V29 13 cos5. (153) 
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Note that the atmospheric sum rule with accurately maximal atmospheric mixing implies 
maximal CP violation with 5 ~ 7r/2 as we observed numerically. On the other hand, 
non-maximal atmospheric mixing is linked to non-maximal CP violation. For 77 = 7r/3 
the predictions shown in Table [5] for the small deviations of the atmospheric angle from 



maximality are well described by the sum rule in Eq. 153 Clearly the atmospheric sum 
rule can be tested in the proposed future high precision neutrino experiments. 
Using the deviation parameters in Eq|2j the above results may be expressed as, 

s pa 0, r pa pa A, S pa ±7l72, a par cos 5, (154) 

3m 3 

corresponding to approximate TBC mixing. Note the novel connection between the 
Wolfenstein parameter and the ratio of neutrino masses, which, if postulated to be an 
exact relation, 

4 mo 
A = -— , 155 

3m 3 

would predict m 2 /m 3 ~ 0.17. 

The TBC mixing matrix in Eqjwwith 5 ~ ±7r/2 becomes, 



1(1 -^ 2 ) U 1 - 1 ^) ^73 A ^ 

^bc-I -^(1±<A) ^(iTfA) ^(1-iA 2 ) 



+ 0(X% (156) 



Ted^) -^( ldh f A ) A( 1 -i A2 )/ 
which is then completely determined in terms of A with no free parameters. 

10 Conclusions 

The type I see-saw mechanism provides a beautiful understanding of the smallness of 
neutrino masses as being due to the heavy right-handed Majorana neutrino masses. 
However, despite its attractive features, the see-saw mechanism provides no understand- 
ing of the observed approximate mixing, characterised here as TBC or TBC2 mixing. 
It also provides no insight into either the ordering (i.e. normal or inverted) or the mass 
scale of the neutrinos (i.e. the mass of the lightest neutrino). Moreover the see-saw 
mechanism is difficult to test experimentally, and typically contains more parameters 
than physical observables. Finally, apparently "unnatural" cancellations can take place 
when constructing the effective neutrino mass matrix from the see-saw parameters. 

One attractive idea which avoids the last problem of "unnatural" cancellations in the 
see-saw mechanism is that the right-handed neutrinos contribute sequentially to the light 
effective neutrino mass matrix with hierarchical strength, leading to the prediction of a 
normal mass hierarchy of physical neutrino masses 7713 3> m 2 ^> m^. This is consistent 
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with recent Planck results which provide no evidence of quasi-degeneracy. The two 
right-handed neutrino model then emerges as the limiting case of a three right-handed 
neutrino model with sequential dominance in the limit that mi — > 0. 

We have considered the minimal type I see-saw model with normal neutrino mass 
hierarchy consisting of two right-handed neutrinos with a zero Yukawa coupling of the 
"dominant" right-handed neutrino to the electron neutrino. The number of observables 
is precisely equal to the number of see-saw parameters in this case. It then becomes 
possible to derive a new and exact "master formula" for such a minimal see-saw model 
which relates see-saw parameters to physical neutrino observables. We hence showed 
that sequential dominance follows automatically and a second texture zero is excluded. 
Using the "master formula" we then explored simple ratios of Yukawa couplings in both 
the dominant and subdominant columns which could explain why the bound on the 
reactor angle #13 < m^jvn^ is saturated by the measurement of #13 ~ 0.15 by Daya Bay 
and RENO. 

The results in FigJ2J show that the ratio of Yukawa couplings \z\\ in the dominant 
column need not be precisely equal to unity, but may vary from \z\\ = 0.5 — 1.5, de- 
pending particularly on the CP violating phase 8. The simple ratio of unity seems to be 
associated with maximal oscillation phase 5 ~ tt/2 for TBC which is an interesting result 
of the analysis. Fixing S = n/2 and \zi\ = 1, the results in FigJ3]for TBC mixing then 
show that simple ratios of subdominant Yukawa couplings \z 2 \ and [^3) and invariant 
see-saw phases 772 and 773 are also possible for particular values of the Majorana phase (3, 
with some examples shown in Table [Tj Analogous results for TBC2 mixing allow some 
new possibilities with 5 = and \zi\ = 0.5 which are displayed in Figjijand Table [2j 

From the point of view of model building we have focussed on the simplest possible 
examples of ratios of Yukawa couplings and see-saw phases, later explaining their origin 
in terms of an underlying A4 family symmetry broken by flavons with simple vacuum 
alignments. We have been particularly interested in examples where the two see-saw 
phases are equal, namely r\ = r] 2 = 773 since it is simpler to justify vacuum alignments 
when the components of a particular flavon have equal phases. In this case then there 
is only a single input see-saw phase 77, which controls both leptogenesis and the two 
PMNS phases 5 and /3, providing a direct link between cosmological CP violation re- 
sponsible for matter-antimatter asymmetry, and the CP violation measured in neutrino 
oscillation experiments. This appears to be a new possibility not discussed before, since 
people hitherto have focussed on the case of this link being provided by having two 
texture zeroes with two right-handed neutrinos. It would also be worth studying flavour 
dependent Ni and N 2 leptogenesis in detail for such examples along the lines of [38] . 

We have discovered two particularly simple possibilities based on the atmospheric 
flavon alignment (0,1,1) and a solar flavon alignment (1,1,3) or (1,3,1), involving a 
relative phase 77 = =F7r/3, which we call CSD3. We emphasise the high degree of predic- 
tivity of CSD3 which involves only two input mass parameters to fix the three neutrino 
masses, three mixing angles, one Dirac CP phase and two Majorana phases (nine physical 
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observables) . We have also seen that the baryon asymmetry is proportional to sin(2r/), 
which is the only see-saw phase for CSD3, providing a link between matter-antimatter 
asymmetry in leptogenesis and low energy CP violation in neutrino oscillation experi- 
ments. We have checked numerically that CSD3 provides a good description of PMNS 
mixing parameters. Fixing the relative phase rj = =F7r/3 by spontaneous CP violation, 
the parameters of PMNS matrix then only depend on the neutrino mass ratio 1712/1713 
as shown in Table [5] for the solar flavon alignment (1,3, 1). We find it remarkable that 
such a simple scheme can lead to such a successful description of all current neutrino 
data, corresponding to approximate TBC mixing to within an accuracy of one degree, 
together with predictions for the phases 5 ~ ±7r/2 and /3 ~ =f70°. 

With CSD3 defined as in Eq|94]we have an example of a minimal predictive see-saw 
model with normal neutrino mass hierarchy. There are many ways to test the model 
by the forthcoming neutrino experiments and it is worth going through them. One way 
to exclude the model is to measure the mixing angles to high precision since according 
to Table [5] these lie in very restricted ranges close (but not quite equal) to the TBC 
ansatz. Note the dependence of the reactor angle on the neutrino mass ratio 1712/1TI3, 
which provides another high precision test of the model. Since the model predicts TMi 
mixing, the deviations are in fact predicted analytically. The solar angle is related to the 



reactor angle by the TMx relation in Eq,149, leading to a value of 612 ~ 34°. The TMi 



atmospheric sum rule in Eq jl53| also leads to an analytic understanding of the smallness 
of the deviations of the atmospheric angle from maximality due to the prediction of 
the approximately maximal oscillation phase 5 ~ ±7r/2. Thus, if the atmospheric angle 
is determined to be non-maximal, then the model will be excluded. It is also worth 
recalling that this is a two right-handed neutrino model with a normal neutrino mass 
hierarchy, and so the lightest neutrino mass is predicted to be zero mi = 0. This 
implies that one of the Majorana phases is irrelevant and may be taken to be zero, 
while the remaining Majorana phase is predicted to be rj 70°, although this is hard 
to measure since the neutrinoless double beta becay observable is predicted to be very 
small, mpp rj |2.6 + 1.15e _l ^ +2 ^|.10 -3 eV rj 2.5. 10 -3 eV. Nevertheless this is significant 
since the observation of neutrinoless double beta decay in the near future would exclude 
this model, as would any cosmological indication of partially degenerate neutrinos. 

We have seen that CSD3 can be realised in an indirect model in which some family 
symmetry is completely broken in the neutrino sector, rather than a direct model where 
some subgroup of the family symmetry is preserved in the neutrino sector. We can 
understand why this is the case by recalling that the neutrino mass matrix predicted by 



CSD3 in Eq. 146 depends on two mass parameters which determine the non-zero neutrino 
masses and these are two of the input parameters of the model. Therefore the Klein 
symmetry of the neutrino mass matrix cannot depend on any of the generators of the 
family symmetry, even though an accidental partial Klein symmetry of TMx mixing is 
predicted. The simplicity of CSD3 as an origin of the reactor angle in the framework of 
indirect models contrasts with the direct approach to the reactor angle which typically 
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involves quite large family symmetry groups. Alternatively one may appeal to charged 
lepton or higher order corrections or else follow the semi-direct approach with only 
half the Klein symmetry preserved. We have proposed an indirect A 4 x Z\° family 
symmetry model which, at lowest order and without charged lepton corrections, explains 
the CSD3 atmospheric flavon alignment (0, 1, 1) and the solar flavon alignment (1, 3, 1) 
via a mechanism based on F-term alignment, with r\ = =Fvr/3 fixed by spontaneous CP 
violation. 

In conclusion, we began in a model independent way by considering a two right- 
handed neutrino model with a dominant texture zero, which has the feature that the 
number of physical observables is equal to the number of see-saw parameters, which we 
related via a master formula. Using the master formula, and the most recent oscillation 
data we then searched for the simplest ratios of Yukawa couplings which can explain 
the data, i.e. the minimal see-saw model with a normal neutrino mass hierarchy which 
can predict the PMNS parameters. Our quest was rewarded by the discovery of CSD3 
which is remarkably simple and economical. With only two input mass parameters, 
CSD3 describes nine physical observables, namely the PMNS mixing matrix plus the 
three neutrino masses. We have constructed an A^ x Z\° model based on CSD3, in 
which the input phase rj = =F7r/3 is determined by spontaneous CP violation and the 
PMNS mixing matrix is completely fixed and only depends on the neutrino mass ratio 
vn^jvno,. The model thus fixes all the mixing angles and phases contained in the PMNS 
matrix and corresponds to approximate Tri-bimaximal-Cabibbo mixing, to an accuracy 
of one degree, with 5 ~ ±7r/2. It is a candidate for the most minimal predictive see-saw 
model consistent with current data and can easily be tested by the forthcoming neutrino 
experiments. In particular the prediction of maximal atmospheric mixing and a normal 
neutrino mass hierarchy will be tested quite soon. 
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